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1  Introduction 


The  increase  use  of  explosive  weaponry  in  military  conflict  and  terrorist  at¬ 
tack  in  recent  years,  has  increased  the  number  of  combat  related  injuries 
sustained  by  soldiers  and  civilian  induced  by  the  blasts.  The  rate  of  ocu¬ 
lar  traumatic  injuries  to  all  injuries  during  the  World  War  II  become  ap¬ 
proximately  6  times  larger  in  compare  to  the  Operation  of  Desert  Storm 
HI  EM,  which  make  the  ocular  injuries  as  a  fourth  most  common  military 
deployment  related  injury  [5].  The  traumatic  injuries  due  to  a  blast  event, 
can  be  induced  by  four  potentially  injurious  mechanisms;  primary  from  the 
blast  overpressure,  secondary  from  propelled  fragments,  tertiary  from  blunt 
impact,  and  quarternary  While  secondary,  tertiary,  and  quarternary 

injury  mechanisms  are  readily  diagnosed  and  managed  within  the  military’s 
casualty  care  system,  specific  injury  mechanisms  unique  to  primary  blast  are 
one  of  the  concern  to  ophthalmologists  [8].  The  severity  of  primary  eye  in¬ 
jury  are  correlated  with  the  high-pressure  shock  front  and  the  subsequent 
wave  of  lower  sub-atmospheric  pressure |9] ,  threshold  overpressure  |10j.  re¬ 
flection  of  the  shock  wave  by  the  orbit  HUHEOE],  and  facial  features  around 
the  eye  [IT].  There  is,  however,  a  dearth  of  clinical  data  that  could  verify 
the  hypotheses  above  and  establish  the  mechanism  of  the  injury.  Measur¬ 
ing  and  assessing  the  influence  of  these  factors  is  difficult  because  survivable 
primary  blast  injuries  are  likely  accompanied  by  injuries  from  fragments  and 
blunt  force  trauma  and  are  thus  more  difficult  to  distinguish  and  enumerate. 
Moreover,  the  severity  of  the  blast  injuries  and  distance  of  the  tertiary  care 
facility  from  the  injury  site  means  that  often  patients  are  unable  to  recount 
the  injury  event,  and  witnesses  are  unavailable. 

Hines-Beard  et  al.  (E|  tested  the  effects  of  blast  overpressure  on  mice, 
by  firing  short  bursts  of  pressurized  air  through  paintball  gun  barrels,  which 
inflicted  closed  eye  injuries  with  features  similar  to  those  seen  in  patients 
with  ocular  blast  trauma.  Alphonse  et  al.  [T6]  studied  the  effect  of  low- 
pressure  blasts  from  fireworks  and  gunpowder  charges  on  human  cadaver 
eyes  but  reported  chances  of  only  minor  corneal  abrasion  due  to  material 
projected  from  the  fireworks  during  explosion.  Risk  of  severe,  physiological 
ocular  damage  was  found  to  be  less  than  0.01%  in  this  study  [16].  Sherwood 
et  al.  P2|  studied  the  anatomic  features  of  primary  blast  ocular  injury  ob¬ 
served  by  using  a  postmortem  porcine  eye  model  exposed  to  various  levels 
of  blast  energy.  Bhardwaj  et  al.  [ITJ  demonstrated  the  effect  of  facial  fea¬ 
tures  on  the  blast  loading  of  the  eye  and  evaluated  the  ocular  biomechanical 


response  to  blast  loading  while  taking  into  account  the  propagation  of  the 
acoustic  wave  inside  the  eye  and  the  associated  ocular  deformation.  Due  to 
the  importance  of  the  facial  features  on  the  profile  of  the  pressure  induced 
by  blast  and  challenges  to  performing  the  experimental  studies  in  full  scale, 
computational  modeling  can  improve  the  understanding  of  the  mechanism 
of  the  primary  injuries.  Computational  models  have  been  used  recently  to 
understand  the  mechanism  of  the  primary  blast  injuries  of  the  eye,  induced 
by  blast  overpressure.  Stitzcl  and  Weaver  [T8j  developed  an  in-house  finite 
element  model  for  blast  loading  on  human  eye,  which  used  the  pressure  on 
the  eye  obtained  as  a  function  of  mass  of  TNT  charge  and  its  distance  from 
the  subject.  Rossi  et  al.  [7]  tested  the  dynamic  response  of  the  eye  to  blast 
overpressure  and  aimed  to  explain  damage  to  the  macula  and  optical  nerve 
head  due  to  the  setting  up  of  a  standing  wave  inside  the  ocular  orbit.  Espos¬ 
ito  et  al.  [jT9j  proposed  a  computational  model  that  presented  the  effects  of 
peak  overpressure  and  blast  duration  on  macular  trauma  and  expressed  the 
possibility  of  severe  damage  of  ocular  tissues  subjected  to  compression  and 
traction  loading. 

To  this  end,  a  computational  model  which  include  the  main  internal  oc¬ 
ular  structures,  spatially  varying  thickness  of  cornea-sclera  shell  and  tissue 
properties,  have  been  utilized  to  help  us  in  understanding  the  mechanism  of 
primary  injuries  and  also  provide  a  a  measurable  estimation  of  risk  of  injuries 
by  using  the  risk  of  injuries  developed  by  Brozoski  et  al.  |20|. 


2  Body 


2.1  Methods 


A  coupled  fluid-structure  interaction  (FSI)  solver  [H]  that  employ  the  fol¬ 
lowing  methods  have  been  used: 

•  Finite  difference  compressible  flow  solver  for  the  propagation  of  the  blast 
wave. 

•  Finite  element  elastodynamic  solver  with  finite  deformation  of  the  human 

eye. 

•  Sharp-interface  immersed  boundary  method  for  fluid-structure  interaction. 
The  details  of  the  above  three  main  components  of  the  modeling  can  be 
found  in  ca.  In  this  study  the  blast  was  modeled  for  the  “base  model” 


which  will  be  described  in  section  (2.1.1)  is  considered  as  a  2-kg  TNT  (Trini- 


Figure  1:  a  Computational  domain  and  initial  conditions  of  the  simulation 
for  the  blast  directly  in  front  of  the  face 


trotoluene)  explosion  in  front  of  the  face,  at  a  distance  of  2.5  m  from  the 
face.  These  parameters  are  obtained  based  on  the  conditions  of  field  tests 
for  blast  exposure,  conducted  by  the  US  Army  Research  program  [21].  The 
initial  conditions  for  the  blast  in  all  simulations  were  identical  to  those  used 
hi  P3j.  The  details  of  the  explosion  models  can  be  found  in  |14j.  A  schematic 
of  the  initial  condition  of  the  blast  is  shown  in  the  figure  [TJ 

2.1.1  FE  models  of  the  head  and  eye 

To  characterize  the  deformation  of  the  globe  resulting  from  the  blast  wave, 
the  FE  model  of  the  head  developed  by  Rajneesh  et  al.  0,  in  which  the  eye¬ 
lid  and  skin  were  neglected  (referred  to  as  “skull  model” )  has  been  used.  The 
anthropometric  parameters  of  a  specific  21-year-old  male,  given  in  Weaver 
et  al.  [22]  is  employed.  The  specific  anthropometric  data  used  for  developing 
the  skull  model  can  be  found  in  Ref.  As  it  is  shown  in  Fig.  [2]  the  model 
includes  all  important  internal  features  of  the  eyes  including  scleral,  limbus, 
cornea,  aqueous  humor,  lens,  ciliary  zonule,  ciliary  muscles,  vitreous  humor, 
retina,  choroid,  lamina  cribrosa,  pre-laminar  neural  tissue  and  the  surround¬ 
ing  orbital/fatty  tissue.  The  friction  between  the  various  ocular  components, 
and  the  skull,  were  not  consider  and  they  fitted  together  seamlessly.  This 
model  was  described  in  detail  in  the  annual  report  for  2013. 

The  corneo-scleral  shell  is  the  outer,  fibrous  layer  of  the  eye  that  protects 
the  various  intraocular  components  and  is  the  main  load  bearing  of  the  eyes. 
The  corneo-scleral  shell,  exhibit  significant  regional  variations  in  thickness, 
and  in  order  to  have  accurate  representation  of  the  behavior  of  the  eye  under 
blast  loading,  these  variations  should  be  taken  into  account.  In  current  study, 
the  sclera  is  constructed  based  on  the  data  obtained  from  high-field  micro 


Figure  2:  Ocular  components  in  exploded  view  with  associated  zone  numbers 
as  described  in  Table  [lj  (Inset)  Assembled  view  of  the  eye  depicting  the 
extraocular  tissues,  sclera,  limbus  and  cornea. 


Table  1:  Zone  numbering  of  different  components 


Component 

Zone  number 

Sclera 

1 

Limbus 

2 

Cornea 

3 

Aqueous 

4 

Lens  (Nucleus,  Cortex  ) 

5 

Ciliary  Zonule 

6 

Ciliary  Muscle 

7 

Vitreous 

8 

Retina 

9 

Choroid 

10 

PLNT 

11 

LC 

12 

Extraocular  tissue 

13 

MRI  measurements  performed  by  Norman  et  al.  [23].  The  corneal  thickness 
is  assumed  to  increase  linearly  from  the  apex  to  the  limbus.  In  the  model 
proposed  by  Stitzel  et  al.  121,  corneal  thickness  was  reported  as  0.52  mm 
at  the  corneal  apex  and  0.66  mm  at  the  limbus.  In  order  to  have  a  linear 
transition  in  thickness  between  tissue,  the  limbus  is  modeled  as  part  of  the 
cornea.  The  retinal  and  choroidal  complexes  were  treated  as  shells  of  uniform 
thickness,  terminating  posteriorly  at  the  ONH.  The  retina  runs  anteriorly  up 
to  the  ora  serrata,  while  the  choroid  extends  up  to  the  ciliary  body.  The 
retina  thickness  considered  as  0.24  mm,  which  was  reported  by  Wagner- 
Schuman  et  al.  [2S|  for  Caucasian  males.  The  choroid  thickness  reported  as 
0.35  mm  by  Ikuno  et  al.  |26j.  based  on  measurements  of  mean  sub-foveal 
choroidal  thickness  for  healthy,  Japanese  individuals,  with  a  mean  age  of  39 
years.  The  lamina  cribrosa  (LC)  is  an  extension  of  the  peripapillary  sclera, 
which  is  a  mesh-like,  multilayered  network  of  collagen  fibers  that  fit  into  the 
scleral  canal  wall,  and  allow  the  passage  of  nerve  fibers  and  a  central  retinal 
artery  into  the  globe.  LC  was  modeled  as  a  disc  from  a  spherical  section  of 
thickness  0.45  mm,  based  on  measurements  for  healthy  eyes,  by  Jonas  et  al. 
ra-  The  sections  of  the  retinal  and  choroidal  shells,  adjacent  to  the  LC, 
were  grouped  as  the  prelaminar  neural  tissue  (PLNT).  The  PLNT  is  similar 
to  retinal  tissue,  consisting  of  nerve  fibers.  Similar  to  the  LC,  the  PLNT  was 
modeled  as  a  spherical  section  disc  of  thickness  0.59  mm  (retinal  thickness 
+  choroidal  thickness). 

2.1.2  Material  properties 

There  is  a  dearth  in  the  literature  describing  the  dynamic  tests  and  me¬ 
chanical  behavior  of  ocular  tissue  subjected  to  blast  impact.  To  best  of  our 
knowledge,  literature  containing  the  mechanical  properties  of  ocular  tissue 
subjected  to  blast  impact  experiments  is  very  limited.  In  the  present  study, 
all  ocular  components  were  treated  as  Neo-Hookean  solids  and  the  best  avail¬ 
able  material  properties,  available  in  literature  were  used.  The  density  of  the 
sclera  was  taken  as  1,400  kg/m 3  (Stitzel  et  al.|24j).  and  the  bulk  modulus 
were  calculated  based  on  the  measured  speed  of  sound  in  the  sclera  [28j,  as 
3570  MPa.  The  dynamic  shear  modulus  of  the  sclera  was  taken  as  10  MPa, 
based  on  the  high  strain  rate  measurements  of  Bisplinghoff  et  al.  [22] .  The 
density  of  the  cornea  is  reported  as  1150  kg/m 3  by  Esposito  et  al.|19j  and  the 
bulk  modules  is  reported  as  3020  MPa  by  Steinert  et  al.  [30j.  The  shear  mod¬ 
ulus  of  cornea  obtained  by  inverse  method  is  reported  by  Esposito  et  al.  [19] 


as  1.66  MPa.  The  properties  of  limbus  considered  to  be  the  same  as  cornea. 
The  density  and  bulk  modulus  of  ciliary  zonule  and  ciliary  muscle  considered 
the  same  as  those  reported  for  the  lens.  The  value  for  the  shear  modulus  are 
calculated  based  on  the  reported  elastic  modulus  by  Van  Alphen  et  al.  |3T| 
as  119  Pa.  The  density  of  the  aqueous  and  vitreous  humor  were  reported 
by  Duck  et  al.  [32]  as  1,003  and  1,009  kg/m3,  respectively.  Shear  modulus 
for  the  vitreous  humor  considered  as  7  Pa,  based  on  the  quasi-steady  mea¬ 
surements  of  the  shear  modulus  of  bovine  vitreous  by  Nickerson  et  al.  [33] 
and  the  shear  modulus  of  aqueous  considered  the  same  as  that  of  vitreous. 
The  density  of  both  nucleus  lens  and  cortex  lens,  considered  as  1100  kg/m3 
(Esposito  et  al.  pH?])  and  the  bulk  modulus  of  both  nucleus  lens  and  cortex 
lens  is  calculated  based  on  the  measured  speed  of  sound  in  the  lens  (Steinert 
et  al.  01  )as  2950.  Shear  modulus  of  nucleus  lens  and  cortex  lens  was  taken 
as  39  Pa  and  98.3  Pa,  from  the  [19].  Retinal  density  is  considered  as  1100 
kg/m3,  as  reported  by  Esposito  et  al.  pj9].  Due  to  lack  of  veritable  data, 
and  since  the  retina  is  composed  mainly  of  cellular  tissues  and  nerve  fibers, 
the  retinal  bulk  modulus  was  approximated  the  same  as  water  (2265  MPa). 
The  Young’s  modulus  of  retinal  tissues  exhibits  large  variations(  PH  135])- 
The  shear  modulus  of  retina  is  considered  as  3.51  kPa  which  is  average  value 
of  the  estimated  shear  modulus  of  porcine  retina  (3.35-3.67  kpa)  with  using 
contact  angle  measurement  of  a  spreading  droplet  reported  by  Grant  et  al. 
[36].  The  material  parameters  for  choroid  considered  as  those  for  the  retina. 
The  Young’s  modulus  for  lamina  cribrosa  was  reported  as  0.3  MPa  by  Nor¬ 
man  et  al.  63.  considering  Poisson’s  Ratio  of  0.47,  the  shear  modulus  was 
calculated  as  0.102  MPa.  Density  of  the  PLNT  considered  the  same  as  the 
density  of  retina.  The  Young’s  modulus  and  Poisson’s  ratio  of  the  PLNT  ob¬ 
tained  based  on  mean  value  of  the  ranges  prescribed  by  Sigal  [38]  which  lead 
to  the  shear  modulus  as  17.3  kPa.  The  properties  for  the  extraocular  tissues 
and  orbital  fat  were  taken  as  those  reported  for  orbital  fat.  The  density  and 
bulk  modulus  were  assumed  to  be  equal  to  those  of  water  1000  kg/m3  and 
2202  MPa,  respectively  ([36]  00]  01]  )•  The  shear  modulus  of  orbital  fat  was 
taken  as  900  kPa  from  the  measurements  of  Schoemaker  et  al.  [42]  for  simian 
retrobulbar  tissues.  The  properties  of  the  ocular  components,  used  in  our 
simulations,  are  as  summarized  in  Table  [2] 


Table  2:  Material  properties  used  for  the  various  ocular  components. 


w  water,  r  retina,  s  sclera,  1  lens 


Component 

p  ( kgm  3) 

Ref 

K  (MPa) 

Ref 

//  (KPa) 

Ref 

Sclera 

1400 

(231 

3571 

1000 

m 

Limbus 

1150 

ra 

3020 

m 

1660 

ra 

Cornea 

1150 

ra 

3020 

ra 

1660 

ra 

Aqueous 

1003 

ra 

2266 

El 

0.007 

0i 

Nucleus  Lens 

1100 

ra 

2950 

0i 

0.039 

na 

Cortex  Lens 

1100 

ra 

2950 

M 

0.0983 

na 

Ciliary  Zonule 

1000 

1 

2950 

1 

0.119 

ra 

Ciliary  Muscle 

1100 

1 

2950 

1 

0.119 

c 

Vitreous 

1009 

m 

2353 

0i 

0.007 

ra 

Retina 

1100 

ra 

2200 

w 

3.51 

Choroid 

1100 

ra 

2200 

w 

3.51 

ra 

PLNT 

1100 

r 

3571 

s 

17.3 

LC 

1400 

s 

3571 

s 

102 

01 

Extraocular  tissue 

1100 

W 

2202 

8 

500 

m 

2.1.3  Results  and  Discussion 


This  section  is  organized  as  follows.  First,  the  results  of  the  simulations  of 
the  skull  model  to  characterize  the  biomechanical  deformation  of  the  the  eye 


caused  by  blast  pressure  loading  are  presented.  In  Section |2.1.4[  the  possible 
injuries  caused  by  blast  over  pressure  which  are  reported  in  literature  are 
correlated  with  the  maximum  stress  obtained  in  different  ocular  components 
and  also  the  risk  of  injuries  are  calculated  based  on  the  injury  risk  available 
in  literature. 

The  maximum  principal  stress  and  maximum  intraocular  pressure  (IOP) 
show  periodic  behavior  (Fig.  [4]  and  Fig.  [3]).  The  overall  increase  in  maxi¬ 
mum  principal  stress  and  IOP  are  due  to  the  time-varying  pressure  loading 
boundary  condition,  and  the  periodic  behavior  of  the  stress  and  IOP  for  oc¬ 
ular  components  forming  the  posterior  wall  of  the  globe  can  be  explained  by 
the  wave  propagation  inside  the  orbit.  Furthermore,  in  order  to  verify  this 
hypothesis,  a  first  order  time  scaling  analysis  of  the  wave  propagation  was 
performed,  the  detail  can  be  found  in  |14j.  The  time  period  calculated  from 
the  simulation  was  around  0.06  ms  (  Fig.  [d]and  Fig.  [3])  which  is  the  same 
order  of  magnitude  as  the  time  period  calculated  by  the  scaling  analysis.  Oc- 


Time  (ms) 


Figure  3:  Time- varying  maximum  intraocular  pressure  (IOP). 


ular  components  forming  the  posterior  wall  of  globe,  experience  higher  oscil¬ 
latory  tensile  and  compressive  stress  (Fig  [3])  and  ocular  components  forming 
the  anterior  structures  of  globe  quickly  reached  their  peak  principal  stress 
and  quickly  decayed(Fig[5]).  Both  trend  of  stress  behavior  agree  with  results 
reported  by  Rossi  et  al.  [7]. 

The  pressure  inside  the  eye  orbit  was  calculated  as  the  mean  of  the  three 
components  of  the  principal  stress.  The  maximum  principal  stress  (Si)  and 
maximum  pressure  in  1.38  ms  are  shown  in  the  transverse  plane  of  the  eye 
orbit  in  Fig.  respectively.  The  highest  maximum  principal  stress  is 

in  thinnest  part  of  sclera  and  the  highest  pressure  can  be  noticed  near  the 
orbital  wall  due  to  the  wave  scattering  by  the  orbital  wall.  The  IOP  was  the 
highest  in  the  posterior  region  of  the  vitreous.  The  contours  of  the  von  Mises 
stresses  in  the  transverse  plane  are  plotted  in  [7]-c,  in  which  the  largest  von 
Mises  stresses  appeared  in  the  thinnest  part  of  the  sclera  and  near  the  optic 
nerve  head.  As  it  can  be  seen  in  Fig.  [7]-c,  the  sclera  is  distorted  on  the  right 
side,  due  to  asymmetric  pressure  loading. 
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Figure  4:  High  amplitude  time  varying  maximum  principal  stress 


Figure  5:  Low  amplitude  time  varying  maximum  principal  stress 
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Figure  6:  von  Mises  stresses  shown  in  the  transverse  and  sagittal  plane  for 
case  3 


Figure  7:  (a)  Maximum  principal  stress  in  the  transverse  plane  of  the  or¬ 
bit,  (b)  Maximum  pressure  in  the  transverse  plane  of  the  orbit,  (c)  Maximum 
von  misses  stress  in  the  transverse  plane  of  the  orbit. 


2.1.4  Implications  to  the  eye  injury 

As  shown  in  Fig.  [3] for  the  baseline  case,  the  maximum  IOP  inside  the  globe 
reached  0.22  MPa,  corresponding  to  around  1650  mmHg,  which  is  two  orders 
of  magnitude  larger  than  the  15  mmHg  physiologic  IOP  for  a  healthy  eye. 
The  von  Mises  stresses  were  calculated  largest  in  the  sclera  wall  at  the  site 
of  muscle  attachments,  which  indicates  possible  risk  of  muscle  tearing  which 
agree  with  recent  experiments  of  blast  exposure  on  rats  which  have  reported 
torn  extraocular  muscles  (Hines-Beard  et  al.|15j).  In  addition,  as  shown  in 
Fig.  [6]  the  maximum  von  Mises  stress  in  sclera,  choroid  and  ciliary  zonule 
show  the  possibility  of  internal  scleral  delamination  (Sherwood  et  ah  [17]). 
chorioretinal  detachments  (Rossi  et  ah  [7],  Sherwood  et  ah  dl ).  and  lens 
dislocation  (Stitzel  et  ah  [24] ) . 

We  assessed  the  risk  of  eye  injury  using  the  functions  published  by  Bro- 
zoski  et  ah  [20]  for  blunt  impact  on  the  eye  by  a  projectile.  The  injury  risk 
function  based  on  survival  analysis  with  maximum  likelihood  is  expressed  as: 


Injury  risk 


100% 


(1) 


where  x  is  projectile  normalized  energy  (kJ/m2),  f3  is  a  dimensionless  param¬ 
eter  that  depends  on  the  injury  type,  and  parameter  a  is  expressed  in  terms 
of  m2 /kJ .  The  values  of  a  and  /3  are  obtained  from  the  maximum  likelihood 
parameter  estimation  technique  for  different  types  of  eye  injuries  which  were 
reported  by  Brozoski  et  ah  ra-  Unfortunately,  the  authors  did  not  report  a 
and  /3  for  corneal  abrasions.  Consequently,  the  injury  risk  function  based  on 
logistic  regression  which  is  also  developed  by  Brozoski  et  ah  [20]  were  used  to 
assess  the  risk  for  corneal  abrasions: 


Injury  risk 


1 

1  +  ea-Px 


100% 


(2) 


The  maximum  IOP  obtained  in  our  simulation  with  the  projectile  normalized 
energy  for  blunt  impact.  The  correlation  between  the  projectile  normalized 
energy  and  IOP  was  reported  by  Duma  et  ah  pH]  for  impact  of  different 
projectiles.  We  chose  the  correlation  for  the  impact  of  11.16  mm  aluminum 
projectile  because  it  is  realistic  representation  of  the  exposed  area  of  the  eye 
that  is  affected  by  blast.  The  calculated  percentage  of  injury  risks  corre¬ 
sponding  to  the  maximum  IOP  value  of  0.3  MPa  for  our  simulation  (baseline 


Tabic  3:  Injury  risks  calculated  for  model  with  uniform  thickness  (Uniform) 
and  model  with  non-uniform  thickness  based  on  maximum  IOP  (MPa)  using 
risk  functions  reported  by  [20],  Cornea  Abrasions  calculated  using  logistic 
regression  technique 


Cases 

Corneal  Abrasions 

Hyphema 

Lens  damage 

Retinal  damage 

Globe  rupture 

Max  IOP 

Uniform 

5 

0.2 

<0.01 

<0.01 

<0.01 

0.22 

Non-Uniform 

4 

0.1 

<0.01 

<0.01 

<0.01 

0.19 

case)  are  approximately  4,  0.1,  0,  0,  and  0%  respectively,  for  corneal  abra¬ 
sions,  hyphema,  lens  dislocation,  retinal  damage,  and  globe  rupture. 

The  injury  risk  values  for  all  cases  are  summarized  in  Table  [3] 


2.1.5  Sensitivity  of  the  model  to  sclera  thickness 

In  this  section,  in  order  to  show  the  sensitivity  of  the  response  of  the  model 
with  respect  to  variational  change  in  thickness  of  the  cornea-sclera  shell,  a 
comparison  between  the  non-uniform  and  model  with  uniform  sclera  thick¬ 
ness  is  performed.  The  model  with  uniform  thickness  have  constant  1  mm 
thickness  for  the  sclera.  The  comparison  of  time-varying  corneo-scleral  max¬ 
imum  principalscleral  maximum  principal  stress  for  non-uniform  model  and 
model  with  uniform  thickness  is  shown  in  Fig.  [9]  The  model  with  non- 
uniform  thickness  showed  8%  increase  in  the  maximum  value  of  maximum 
principal  stress  and  3.5%  increase  in  the  average  of  maximum  principal  stress. 
The  comparison  of  time-varying  maximum  IOP  for  non-uniform  model  and 
model  with  uniform  thickness  is  shown  in  Fig.  [8]  The  model  with  non- 
uniform  thickness  showed,  14%  decrease  in  the  maximum  value  of  maximum 
IOP  and  3.3%  decrease  in  the  average  of  maximum  IOP  in  compare  to  uni¬ 
form  model.  The  injury  risk  values  for  the  uniform  model  are  summarized  in 
Table  [3}  Even  though,  the  maximum  cornea-sclera  shell  stress  and  maximum 
IOP  in  both  model  have  very  similar  trends  in  both  models,  however,  as  it 


is  shown  in  Fig.  10  (a),  the  average  of  maximum  von  Mises  stress  in  cornea- 
sclera  shell  for  non-uniform  model  is  higher  (25%)  than  uniform  model.  In 
addition, the  average  of  maximum  von  Mises  stress  in  ciliary  and  choroid  for 
non-uniform  model  are  higher  6%  and  40  %,  respectively  and  4%  lower  in 
limbus. 


Figure  8:  Comparison  among  four  cases  with  different  thickness  for  time 
varying  maximum  intraocular  pressure,  IOP 


Figure  9:  Comparison  among  four  cases  with  different  thickness  for  time 
varying  and  corneo-scleral  principal  stress,  si 


Figure  10:  Comparison  of  time-varying  maximum  von  Mises  stress  for 
(a)sclera,  (b)  limbus,  (c)  ciliary  body,  (d)  choroid 

3  Limitations  of  the  computational  model 

In  the  present  computational  model,  the  tissues  of  the  cornea  and  sclera 
were  assumed  to  be  isotropic  and  spatially  homogeneous  in  the  Neo-Hookean 
constitutive  model,  although  Aghamohammadzadeh  et  al.  [53]  have  demon¬ 
strated  the  anisotropic  and  non-linear  behavior  of  the  sclera.  We  are  cur¬ 
rently  working  on  adding  the  anisotropy  to  the  model.  Although  we  have 
used  the  best  available  material  properties  for  the  ocular  components,  there 
may  be  uncertainty  in  the  material  properties  at  high  blast  loading  rates. 
We  also  did  not  account  for  the  viscoelastic  nature  of  ocular  tissue  which  can 
influence  the  stresses  and  IOP. 


4  Work  in  Progress 

We  are  currently  working  on  incorporating  collagen  derived  anisotropy  of  the 
cornea  and  sclera  in  the  computational  model  to  investigate  their  effects  on 
the  risk  of  injuries,  deformation,  stress  response  of  the  globe  to  the  blast 


wave. 


4.1  Anisotropic  mechanical  properties  of  the  cornea 
and  sclera 

The  mechanical  properties  of  the  sclera  obtain  from  its  fiber-reinforced  mi¬ 
crostructure,  which  consists  of  a  stack  of  collagen  lamellae,  composed  mainly 
of  unidirectional  type  1  collagen  fibers  embedded  in  a  proteoglycan  rich  ma¬ 
trix  ([46]).  The  collagen  fiber  structure  of  the  posterior  human  sclera  was 
recently  investigated  using  wide-angle  X-ray  scattering  (WAXS)  (Pijanka  et 
al.  [37]).  In  WAXS,  X-rays  passing  through  the  ordered  fibrous  collagen 
structure  are  scattered  to  produce  a  diffraction  pattern.  The  quantification 
of  corneal  and  scleral  collagen  fiber  orientation  from  WAXS  patterns  can  be 
found  in  detail  in  [38].  The  WAXS  measurements  showed  the  existence  of 
a  ring  of  highly  aligned  collagen  in  the  region  closest  to  the  ONH.  the  de¬ 
gree  of  anisotropy  in  the  peripapillary  ring  varied  markedly  (between  2-fold 
and  4-fold)  with  circumferential  position  around  the  scleral  canal  and  weak¬ 
est  circumferential  alignment  is  observed  in  the  superior-nasal  quadrant  m ■ 
The  collagen  structure  in  the  midposterior  sclera  exhibited  large  variations 
in  the  preferred  fiber  orientations  as  well  as  a  low  degree  of  fiber  alignment 
for  this  reason  the  rest  of  the  sclera  was  considered  as  isotropic  material  in 


this  study.  In  our  study,  as  it  is  shown  in  Fig.  [XT]  the  peripapillary  sclera  is 
divided  to  4  quadrant  and  different  degree  of  anisotropy  is  assigned  to  each 
quadrant. 

The  material  parameters  for  the  cornea  (/i,  (3,  a  )  is  selected  in  way  to 
match  the  uniaxial  tests  performed  by  Zeng  et  al.  on  human  cornea  [45], 
The  arrangement  of  collagen  fibrils  is  important  to  determine  the  mechan¬ 
ical  strength  of  the  cornea.  X-ray  scatter  intensity  distribution  data  has 
indicated  two  preferred  directions  of  collagen  orientation  at  the  center  of  the 
cornea:  along  the  nasal-temporal  and  superior-inferior  meridians ([50,  I5T]). 
Closer  to  the  limbus,  the  fibers  tend  to  run  in  a  circumferential  direction 
(sa  In  the  center  of  the  cornea,  preferential  orientations  are  dominant  along 
the  superior-inferior  SI  and  the  nasal-temporal  NT  meridians.  At  the  limbus, 
where  the  cornea  and  the  white  sclera  merge,  most  collagen  fibrils  are  aligned 
circumferentially  clearly  forming  a  reinforcing  annulus.  The  anisotropic  ma¬ 
terial  parameters  are  summarized  in  Table  [4] 


Figure  11:  Cornea-sclera  shell  sections 


Table  4:  Anisotropic  material  parameters  of  cornea-sclera  shell,  MA,  mid¬ 
posterior  and  anterior  sclera;  PP,  peripapillary  sclera;  SN,  superior-nasal 
quadrant;  ST,  superior-temporal  quadrant;  IN,  inferior-nasal  quadrant;  IT, 


inferior-temporal  quadrant. 


Component 

/j,(KPa) 

a(Pa) 

P 

Ref 

Fiber  Anisotropy 

Fiber  Orientation 

Ref 

Sclera-MA 

160 

5.00E+02 

215 

11531 

0 

0 

sisij 

Sclera-PP-SN 

160 

5.00E+02 

215 

m 

0.45 

0.74 

GS3J 

Sclera-PP-ST 

160 

5.00E+02 

215 

f52j 

0.68 

0.74 

m 

Sclera-PP-IN 

160 

5.00E+02 

215 

II 5211 

0.59 

0.74 

m\ 

Sclera-PP-IT 

160 

5.00E+02 

215 

II 5211 

0.54 

0.74 

m\ 

Limbus 

60 

1.04E+05 

4.45 

Estimated 

0.37 

0 

m3] 

Cornea 

60 

1.04E+05 

4.45 

Estimated 

0 

0 

TEST 

4.2  Anisotropic  hyperelastic  constitutive  model 

The  sclera  and  cornea,  limbus  was  idealized  as  a  continuum  composite  con¬ 
sisting  of  fiber  families,  each  representing  a  collagen  lamella  with  an  in-plane 
orientation  angle,  0,  embedded  in  an  isotropic  matrix.  The  collagen  structure 
was  described  by  a  continuous  probability  density  distribution  of  the  fiber 
orientations,  D  (</>).  The  matrix  and  collagen  fibers  deform  with  the  tissue  ac¬ 
cording  to  the  macroscopic  deformation  gradient  F,  which  maps  a  vector  from 
the  undeformed  configuration  to  the  deformed  configuration.  The  invariants 
of  the  right  Cauchy  green  tensor  (C  =  F7F  )  are  defined  as  follows: 


h  =  tr(C)  (3) 

h  =  i(tr(C)2  -  tr(C2)  (4) 

h  =  det(C)  (5) 

The  strain  energy  for  component  with  this  material  model  can  be  decomposed 
into  a  matrix  term  and  an  anisotropic  term  dependent  on  the  fiber  stretch 
which  can  be  described  as  follow: 

w(  C)  =  Wmatrix(h(C),I2(C),I3(C))  +  Waniso{  C)  (6) 

The  collagen  fibers  assumed  to  be  identical  and  can  be  described  by  using 
the  same  strain  energy  density 

Wfiber(\f)  =  ^[exp(/3(Aj  -  1))  -  p\}\  (7) 

Where  4 a/3  is  stiffness  of  the  fiber  family  and  (3  is  a  strain  stiffening  parameter 
and  Xf  is  the  stretch  of  the  fiber  determined  by  C  :  N(0),  where  N  is  the 
unit  vector  describing  the  fiber  orientation.  The  anisotropic 

Waniso( C)  =  r  Wfaer(Xf(<j>))D{<j>)d<f>  (8) 

J  —TV 

The  matrix  phase  was  modeled  as  an  isotropic, quasi-incompressible  Neo- 
Hookean  material. 

Wmatrix(h(C),  h (CO)  =  |(/1  -  3)  +  |(/3  -  In (/3)  -  1) 

where  /i  is  shear  modulus  and  k  is  bulk  modulus. 


(9) 


Table  5:  The  range  of  shear  modulus  for  isotropic  parameters  sensitivity 


analysis 


Component 

Upper  range  (Pa) 

Lower  range(Pa) 

Sclera 

106 

107 

Cornea 

1660xl03 

1660  xlO4 

Nucleus 

25.2 

52.8 

Cortex 

33.8 

162.8 

Retina  and  Choroid 

2800 

4200 

Ciliary  Zonule 

95.44 

143.16 

4.2.1  Sensitivity  of  the  model  to  the  material  properties  and  blast 
conditions 

Sensitivity  analysis  will  be  performed  to  investigate  the  effect  of  material 
parameters  and  blast  conditions  (mass  and  distance  of  explosive)  on  stress 
response,  maximum  IOP  and  risk  of  injuries.  The  sensitivity  analysis  will  be 
performed  as  follow. 

Isotropic  We  are  conducting  a  parameter  study  to  understand  the  effect 
of  material  properties  of  the  sclera,  cornea,  lens  (nucleus,  cortex),  ciliary 
zonule,  retina  and  choroid  for  isotropic  model.  The  bulk  moduli  and  density 
were  not  varied.  The  shear  moduli  of  these  components  are  considered  as 
a  variable  in  our  sensitivity  analysis.  The  shear  modules  of  nucleus,  cortex 
and  the  shear  modules  of  retina,  choroid  considered  to  change  together  which 
reduce  the  number  of  variables  of  this  study  to  5  variables.  In  order  to 
understand  the  effect  of  shear  moduli  on  the  maximum  IOP  and  risk  of 
injuries,  a  fractional  factorial  design  is  used  to  construct  a  23  experiments 
which  achieves  resolution  III  design.  The  range  of  the  variation,  for  nucleus 
and  cortex  considered  as  [25.2  52.8]Pa  and  [33.8  162.8]Pa,  respectively,  based 
on  the  values  reported  by  [19].  The  shear  modulus  of  the  corneo  and  sclera, 
varied  from  their  base  value  to  the  10  fold  higher  which  is  in  the  order  of 
the  scleral  modulus  measured  for  high-pressure  rates  [29].  The  variation  in 
shear  modulus  of  ciliary  zonule,  retina  and  choroid  considered  as  20%  of  their 
nominal  values  which  was  reported  in  Table  [2}  The  range  of  shear  modulus 
are  summarized  in  Table  [5] 


Table  6:  Range  of  anisotropic  material  parameters  of  cornea  sclera  shell,  MA, 
midposterior  and  anterior  sclera;  PP,  peripapillary  sclera;  SN,  superior-nasal 
quadrant;  ST,  superior-temporal  quadrant;  IN,  inferior-nasal  quadrant;  IT, 
infer ior-temporal  quadrant . _ _ _ 


Component 

P 

a(3  (KPa) 

Fiber  Anisotropy 

Sclera-MA 

Sclera-PP-SN 

Sclera-PP-ST 

Sclera-PP-IN 

Sclera-PP-IT 

[172  258] 

[86  129] 

0 

[0.37  0.52] 

[0.56  0.82] 

[0.49  0.70] 

[  0.45  0.64] 

Limbus 

Cornea 

[3.56  5.34] 

[370  555] 

[0.33  0.40] 

0 

Anisotropic  A  parameter  study  will  be  performed  to  investigate  the  effect 
of  anisotropic  material  parameters  of  the  cornea-sclera  shell  on  the  stress 
response  and  estimated  risk  of  injuries.  The  bulk  moduli  and  density  were 
not  varied.  As  it  was  described  in  section  |4.1[  different  anisotropic  param¬ 
eters  (preferred  fiber  orientation,  degree  of  anisotropy)  have  been  assigned 
to  different  section  of  cornea-sclera  shell  based  on  WAXS  data.  In  order 
to  understand  the  effect  of  anisotropic  material  parameters  of  cornea-sclera 
shell  on  the  stress  response  and  risk  of  injuries,  a  set  of  simulations,  designed 
based  on  resolution  III  fractional  factorial  analysis,  will  be  performed.  The 
range  of  anisotropic  material  model  for  cornea  sclera  shell  are  summarized 
in  Table  |6l 

Blast  conditions  The  effect  of  the  parameters  of  the  blast  on  the  injury 
risk  and  stress  response  in  ocular  component  will  be  investigated  as  part  of 
sensitivity  analysis.  The  simulation  will  be  performed  with  0.64  and  6.4kg 
TNT  blasts  and  the  result  will  be  compared  to  the  base  model.  Based  on 
the  weight  of  the  explosive  charge  used,  the  standoff  distance  (radius  of  the 
circle)  was  varied,  such  that  conditions  for  0%  lethality  and  pulmonary  injury 
threshold  were  emulated.  [S3] 


4.2.2  Risk  of  Injuries 


We  are  currently  applying  experimental  data  for  blunt  impact  to  asses  the 
injury  risks  associated  with  the  blast  pressure.  This  likely  leads  to  erroneous 
injury  risk  criteria  because  the  IOP  and  stresses  generated  in  the  intraocular 
tissues  by  the  propagating  blast  wave  is  very  different  for  the  same  impact 
energy  comparing  blast  to  blunt  impact.  Blast  events  exposes  the  tissues  to 
significantly  higher  rates  and  the  high  rate  response  of  ocular  tissues  has  been 
shown  to  be  different  than  for  low  rates.  There  is  currently  no  experimental 
data  published  for  primary  blast  injury  to  ocular  tissues.  Sherwood  and 
coauthors  El  HZ!  have  begun  to  examine  damage  in  cadaveric  porcine  eyes 
exposed  to  blast  pressures  generated  by  a  shock  tube.  In  future  work,  we  will 
used  their  data  to  estimate  the  probability  of  sub-globe-rupture  injuries  under 
various  blast  pressure-impulse.  Blast  parameters  will  be  used  to  estimate  the 
likelihood  of  a  specific  type  of  injuries.  A  score  based  system  will  be  used  to 
determine  the  levels  of  trauma  based  on  the  blast  pressure  profile  int 

5  Key  Research  Outcomes 

1.  Developed  a  more  detailed  model  of  the  human  eye  with  internal  ocular 
structures  and  spatially  varying  thickness. 

2.  IOP,  VMS  Mises  stress  and  principal  stresses  were  obtained  in  different 
component  of  the  eye. 

3.  Risk  of  injuries  were  used  to  obtain  the  likelihood  of  specific  injuries 
due  to  the  blast. 

4.  The  findings  of  our  simulations  support  experimental  observations  of 
injury  to  choroid,  sclera,  ciliary  zonule  and  ONH. 

5.  The  effect  of  local  thickness  variation  have  been  investigated  on  the 
stress  response,  internal  ocular  pressure  and  estimated  risk  of  injuries. 

6.  Began  to  include  the  anisotropic  material  model  to  the  cornea-sclera 
shell  in  order  to  describe  the  nonlinear,  anisotropic  stress  response  of 
the  load  bearing  components  of  eye. 

7.  Began  to  estimate  the  effect  of  material  parameters  and  blast  conditions 
in  the  stress  response  of  the  eye  and  estimated  risk  of  injuries 


6  Conclusion 


A  FSI  computational  model  was  presented  to  investigate  the  effect  of  the 
blast  wave  on  the  human  eye.  The  transient  flow  fields  and  pressure  loading 
profile  on  the  human  eye  were  obtained  by  simulating  the  three  dimensional 
blast  wave  propagation.  The  possible  eye  injuries  and  associated  injury  risks 
were  estimated  by  using  the  calculated  intra-ocular  pressure  (IOP)  and  von 
Mises  stress  in  the  various  ocular  components.  The  sensitivity  of  the  model 
to  local  variation  of  the  sclera  were  determined  and  discussed. 

The  future  goals  are  to  consider  anisotropic  material  of  the  cornea-sclera 
shell  and  determine  the  effect  of  anisotropy  on  the  stress  response  of  the  eye 
and  estimated  risk  of  injuries  in  our  model.  In  addition,  we  will  perform  the 
parameters  study  (material  parameters,  blast  conditions)  to  determine  the 
effect  of  each  parameters  on  the  stress  and  estimated  risk  of  injuries.  A  new 
damage  criteria  based  on  the  blast  parameters  will  be  used  to  estimate  the 
likelihood  of  a  sub-globe-rupture  injuries. 
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